Current percolation and anisotropy in polycrystalline MgB 2 
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The influence of anisotropy on the transport current in MgE>2 polycrystalline bulk samples and 
wires is discussed. A model for the critical current density is proposed, which is based on anisotropic 
London theory, grain boundary pinning and percolation theory. The calculated currents agree 
convincingly with experimental data and the fit parameters, especially the anisotropy, obtained 
from percolation theory agree with experiment or theoretical predictions. 
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Soon after the discovery of superconductivity in MgB 2 
0, it was found that grain boundaries do not limit the 
current flow in polycrystalline samples Thus, sin- 

tered bulk samples and wires could be prepared with 
critical current densities appropriate for technical appli- 
cations at low magnetic fields |4|. With increasing field 
the critical current rapidly decreases, becoming zero at 
a field far below the upper critical field. Anisotropy, ob- 
served on thin films |j| for the first time, was identified as 
a possible reason for this strong field dependence . The 
aim of this letter is to evaluate the influence of anisotropy 
on the field dependence of J c in terms of percolation the- 
ory 0, @- Percolation theory had been employed frc- 
uently to analyze currents i n g ranular superconductors 
"ISlIlHIIiQHHOj, especially in high tem- 
perature superconductors, where grain boundaries limit 
the current flow, and also to interpret the irreversibility 
line in MgB2 coatings . We present a percolation model 
for the limitation of the transport current by anisotropy. 
While some details refer to MgB2, the model is applica- 
ble to any granular anisotropic superconductor or, more 
generally, to other systems with spatially inhomogenous 
transport properties. 

The superconductor is assumed to consist of randomly 
oriented grains with exactly identical properties. Ther- 
mal activation of the flux line lattice is neglected, since 
its influence in granular MgB2 should not be much larger 
than in the low-T c materials NbTi or NbaSn @ (this sim- 
plification might be the reason for some systematic devia- 
tions at high temperatures, cf. below) . At fixed tempera- 
ture the upper critical field only depends on the angle be- 
tween the boron layers in the grain and the applied mag- 
netic field according to the anisotropic Ginzburg Landau 
relation (Isf 



B c2 {9) 



B c2 (ir/2) 



sff cos 2 ((9) + sin 2 (<2) 



(1) 



7 denotes the anisotropy factor of the upper critical field, 
i.e. B c2 (ir/2)/B c2 (0) =: b\ 2 /B^ 2 , and 9 is the angle be- 
tween the applied field and the c-axis. This relation was 
confirmed experimentally in MgB2 by torque measure- 
ments 0, We start by considering the resistive 



transition in decreasing temperature at fixed field Bq. 
At a certain temperature, b\ 2 becomes equal to the ap- 
plied field Bq and the resistivity starts to decrease. This 
point is per definition B c2 of the whole sample. When 
the temperature decreases further, an increasing number 
of grains becomes superconducting, depending on their 
orientation. During this transition the system can be 
considered as a mixture of normal and superconducting 
grains. Systems consisting of two materials with differ- 
ent resistances (conducting /insulating 0, 112, El, 0] or 
normal/superconducting .25, 26]) have been investigated 
experimentally and are well understood in the framework 
of percolation theory. In mixtures of normal and super- 
conducting powders the resistivity remains finite as long 
as the probability p, that a grain is superconducting (i.e. 
the fraction of superconducting grains), is smaller than 
the critical probability or percolation threshold p c . The 
percolation threshold strongly depends on the number 
of connections to neighboring grains and, therefore, on 
the preparation conditions. When the fraction of super- 
conducting particles becomes equal to the critical proba- 
bility, the first continuous superconducting current path 
occurs and the resistivity disappears for sufficiently small 
currents. Since the upper critical field increases mono- 
tonically with 9, an angle 9 C corresponds to the critical 
probability within the transition region. If the boron 
planes are randomly oriented, the directions of the c- 
axes are equally distributed in 3D space and their angles 
to any fixed direction are distributed as s'm(9). Integra- 
tion leads to 9 C = arccos(p c ). With B c2 (ir/2) = Bo and 
B c2 (9 c ) = Bq for the onset and the offset of the transi- 
tion, and assuming ^§fr- and 7 to be constant within the 
transition region, the transition width is obtained from 
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AT a represents the broadening of the transition only due 
to the anisotropy and is, therefore, additive to other rea- 
sons for a finite transition width, e.g. material inhomo- 
geneities, thermal fluctuations, thermal gradients in the 
sample, surface effects etc. The numerator of ©, ftw, 
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FIG. 1: Onset and offset (a) of the resistive transition in bulk 
MgB2 and normalized transition width (b). 



can be expressed explicitly and interpreted as a normal- 
ized transition width. It is a monotonically increasing 
function of 7 and p c . 

The resistive transition of a sintered bulk sample Q 
was measured at fixed magnetic field at a current den- 
sity of 2 kA/m 2 . The onset of the transition was defined 
at 95 % of the normal state resistance, the offset at 5 %. 
The results are plotted in Fig. From the difference of 
onset and offset temperature, f tw is calculated and plot- 
ted in Fig. ^3 (black squares) as a function of the mean 
value of these two temperatures. Since the contribution 
of anisotropy tends to zero at zero field, it cannot be 
dominant at low fields, which explains the upturn of the 
normalized transition width at high temperatures. As a 
simple correction, we subtract the zero field transition 
width (0.25 K) at all fields, which leads to the solid line 
(without symbols) in Fig. \\]p. The upturn of the original 
data is removed and f tw becomes nearly temperature in- 
dependent below 30 K, but decreases rapidly at higher 
temperatures. The only temperature dependent factor 
in f tw is the anisotropy and, indeed, a rapid decrease of 
anisotropy near T c was found [2fJ, |27| • However, the as- 
sumption that all other influences on the transition width 
are temperature and field independent, are not appropri- 
ate to draw exact conclusions on the temperature de- 
pendence of the anisotropy. Nevertheless, the increasing 
difference between the onset {B C 2) and the offset (irre- 
versibility field) of the resistive transition with increasing 
magnetic field can be explained by anisotropy and only 
minor corrections are needed to obtain the field depen- 
dence of the transition width observed in experiments. 

In order to derive a model for the critical currents in 
polycrystalline samples, grain boundaries are assumed to 
be the dominant pinning centers. The dependence of 
the critical current density on the ap plie d magnetic field 
(Jc-LB) is, therefore, expected to be |2j| 
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for B < B C 2(0) and zero otherwise. This equation was 
originally developed for isotropic superconductors, but 



anisotropy can be introduced through the angular depen- 
dence of the reversible parameters, k = Ki = B C 2/y2B c 
was used. For the sake of simplicity, we neglect any in- 
fluence of the quality of a specific grain boundary, the 
details of grain boundary pinning as well as the angle 
between the local current and the applied field, which 
may differ due to percolation. The main feature we 
wish to address, is the disappearance of the critical cur- 
rent as (1 — B I B C 2(6)) 2 at high fields and its rapid de- 
crease as B~ 5 at low fields. With these critical cur- 
rent densities of individual grains, we employ results on 
random resistor networks to derive the critical current 
density of the whole sample. There, the current in a net- 
work consisting of equal resistors with a current voltage 
law U oc I a is given near the percolation threshold by 
I = (JeffU 1 / , with the conductivity a e ff oc (p — p c )* 
[29| and t = (d — \)v + (£ — v)/a [3(|, where d denotes 
the dimensionality of the system, £ is about one and the 
correlation exponent v is 0.88 in 3D systems jU- Most 
experimental data were obtained on systems consisting 
of normal conducting (a = 1) and insulating particles 
EH E3- El- In this case, the above power law for 
a e ff was found to be an excellent approximation and t 
to be in reasonable agreement with the theoretical pre- 
diction (t — 1.88). The above dependence of t on a 
was confirmed numerically for 2D nonlinear resistor net- 
works 29]. Its asymptotic (a — ► 00) value, t = (d — l)v, 
was proven for all dimensions (t — 1.76 for d = 3). 
In the following we assume a to be infinite and t to be 
1.76, which makes the current independent of the volt- 
age ("hard" superconductor), but the argument remains 
essentially the same for any other a at fixed voltage. 
Note that the influence of a on t is very weak in 3D 
systems. In a granular superconductor the conductivity 
cr e ff becomes a function of the (local) current density 
J, if p is interpreted as the fraction of grains with crit- 
ical current densities larger than J. With the sample 
cross section tro and with an appropriate normalization 
to get correct results for p = 1, the conductivity (cur- 
rent) becomes I — a e ff — g ( p ^~ Pc ) t J =: a p J. We 



introduce a p , which can be interpreted as the total cross 
section of all current paths existing at a certain fraction 
p of conducting particles. If the current is smaller than 
00 mine( Jc(B, 8)), p is one and the current flows homo- 
geneously through the whole superconductor. Upon en- 
hancing the current by A/, p becomes smaller than one 
and the current density in the remaining current paths 
increases by A J = AI/a p . The remaining capability of 
the grains to carry additional currents is reduced by A J. 
The current can be further enhanced until p = p c . In- 
tegrating over all infinitesimal small AI (or more easily 
A J) and division through ctq leads to the macroscopic 
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critical current density of the sample: 



MB) 



'(B) 
( 



P(J) -Pc 
1 ~Pc 



l dJ 



(4) 



where J™ ax {B) is defined by p(J™ ax {B)) = p c . In or- 
der to evaluate this expression numerically, 1001 grains 
were chosen with sin(8) equally distributed between 
and 1. Their critical currents were calculated according 
to J3J and the probability distribution p was calculated 
basically by counting the number of grains with a crit- 
ical current density above a certain value. The integral 
was then replaced by a summation over the discrete crit- 
ical current densities of the chosen grains. The model 
contains four parameters: the upper critical field B C 2, 
the absolute value of the critical current densities in J2J| , 
the critical probability p c , and the anisotropy 7. While 
B C 2 can be measured directly at high temperatures and 
extrapolated to low temperatures, the remaining param- 
eters are obtained by fitting the logarithm of J C (B) ob- 
tained from I0J to the corresponding experimental data 
with the least squares method. Two samples were used 
to compare experimental results with the predictions of 
the model, a sintered bulk sample and a Cu-sheathed 
powder-in-tube wire, prepared by an in-situ process [33j |. 
The samples were neutron irradiated in a fission reactor, 
in order to change the reversible properties. Details of 
the irradiation experiments can be found in |34i l35| . The 
upper critical field of a similar bulk sample (from the 
same pellet) was determined by the onset of the resistive 
transition (Fig. ^i) and the linear behavior below 30 K 
was extrapolated to lower temperatures. For the wire, 
the onset of the resistive transition cannot be evaluated 
for B c2 due to the low resistance of the copper sheath. 
Therefore, B C 2 was obtained from SQUID measurements. 
The critical current densities of the bulk sample were ob- 
tained from ac susceptibility measurements, at high crit- 
ical current densities by a technique described in [36j |. 
and near the irreversibility line by evaluating the magni- 
tude of the out-of-phase component in terms of the Bean 
model. Since the first method is only applicable, if the 
ac penetration depth is much smaller than the sample 
radius, and the latter, if it is larger, intermediate current 
densities were not evaluated. In the case of the wire, di- 
rect transport measurements were made in liquid helium. 

The results are plotted in Fig. [5^ together with the 
fitted curves (solid lines). Since the excellent agreement 
between experimental and calculated data (obtained by a 
three parameter fit) might not be unique to the proposed 
model, it is most important that the fit parameters (Ta- 
ble [J are indeed in agreement with expectations. The 
anisotropy is found to be 4.4 for both unirradiated sam- 
ples. Literature data range from slightly above 1 to 9 
and are obviously sample dependent. A direct determi- 
nation of the B C 2 anisotropy was made on single crystals 
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FIG. 2: Critical current densities of various samples at low 
temperatures (a) and of bulk MgE>2 at higher temperatures 
(b). The solid lines represent the fit. 



(7 = 4.6 |23), aligned powders (7 = 1.7 37]) and thin 
films (7 = 1.2 — 2 [5|). For sintered polycrystalline sam- 
ples (wires, bulk), the anisotropy was extracted by a spe- 
cial technique from magnetization measurements and es- 
timated to be 6-7 . The large differences indicate that 
disorder affects anisotropy. This is directly confirmed 
here by radiation-induced disorder (Table P). The perco- 
lation threshold p c is similar in the unirradiated samples 
(0.21 and 0.26). Theoretical predictions for site percola- 
tion are about 0.2 for the face centered cubic lattice (co- 
ordination number Z=12) and 0.31 for the simple cubic 
lattice (Z=6) [3!| . Although the grains in our samples are 
not arranged in a regular lattice, our result indicates that 
on average each grain is connected to 6-12 neighboring 
grains. This has been confirmed by SEM investigations. 
Other experimental data range from 0.17 in pressed mix- 
tures of conducting spheres and Teflon powder [24]], to 
0.31 in systems consisting of silver coated and uncoated 
glass beads or 0.47 in granular W-A1 2 3 films 0. 

Experimental data on the unirradiated bulk sample at 
higher temperatures and the fitted J C (B) curves are plot- 
ted in Fig. Eb, the resulting values of 7 and p c in Fig. 
While a temperature dependence of the B C 2 anisotropy 
was also found in single crystals 0, |2(J and powders of 
MgB 2 27], the increase of p c with increasing temperature 
is quite unexpected, since p c should only depend on the 
arrangement of the grains within the sample. In order to 
explain this temperature dependence, we investigate the 
influence of 7 and p c on the critical currents. The dotted 
line in Fig. [2^ represents the calculated current densi- 
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FIG. 3: Temperature dependence of the anisotropy and of the 
percolation threshold in bulk MgB2. 



TABLE I: Experimental data on the upper critical fields 
(Sj 2 ) as weu as anisotropy and percolation threshold obtained 
from the fits of the critical current densities. 



B^ 2 (extrapol.) 
anisotropy 

Bc2 (:= Sj 



bulk (5 K) bulk irr. wire (4.2 K) wire irr. 



h) 



21.4 T 

4.4 
4.9 T 
0.21 



30 T 
2.8 
10.7 T 
0.4 



15.4 T 

4.4 
3.5 T 
0.26 



18.6 T 

2.5 
7.4 T 

0.29 



ties, if p c of the unirradiated bulk sample is enhanced 
by 20 %, leaving the other parameters constant. This 
has very little effect for fields below B^r 2 - An increase 
of the anisotropy by only 10 % (dashed line in Fig. 
significantly decreases the critical currents to much lower 
fields. The very small influence of p c compared to 7 at 
low fields makes the evaluation of the anisotropy quite 
robust against changes in p c . On the other hand, at a 
given anisotropy, p c determines the irreversibility field 
and, therefore, the transition width. The fit parameter 
p c cannot be interpreted only as the percolation threshold 
anymore, because it also contains any other broadening 
effects on the transition. Especially thermal fluctuations 
lead to an artificial increase of the percolation threshold 
at high temperatures, but also material inhomogencitics 
enhance p c . The latter explains the increase of p c af- 
ter irradiation (Table due to the large inhomogeneity 
of the defect structure following neutron irradiation |34| , 
which leads to a significant broadening of the transition 
even at zero field (AT C ). The wire was irradiated in a 
cadmium cover 35] , which reduces the inhomogeneity of 
the resulting defect structure and the enhancement of p c . 

Finally, the normalized transition width ft w was cal- 
culated from the fit parameters of the unirradiated bulk 
sample (Fig. open circles) according to J2J- The 
agreement with the data directly calculated from the 
measured transition width is very satisfactory. We con- 
clude that the width of the resistive transition can be 
quantitatively explained under realistic assumptions for 
the anisotropy (3.1-4.4) and for the percolation proba- 
bility (0.21-0.29), and that this contribution is dominant 
at low temperatures. 

In summary, we have shown that the dependence of the 
critical currents in polycrystalline MgB2 on the applied 
magnetic field and the broadening of the resistive transi- 
tion with increasing magnetic fields can be explained by 
anisotropy and current percolation. The rather low up- 
per critical field B^ 2 , turns out to be responsible for the 
strong field dependence of the critical current densities 
in polycrystalline MgB2. 
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